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abstract 

This paper is the second of a series dealing with the motion of flood waves and 
other waves of translation in open channels. The first paper considered waves 
controlled solely by inertia forces; the present one is an analysis of the combined 
effects of turbulent friction and inertia. The basic equation of motion for gradu- 
ally varied unsteady flow in prismatic channels is derived from fundamental 
principles. The effect of the velocity distribution in the original undisturbed 
current on the motion of short waves is investigated, and the effects of wave 
height, curvature of profile, and fluid friction on the celerity of a wave-volume 
element is analyzed in detail. The deformation of a straight sloping front and the 
change of height of an abrupt wave front is treated. Special emphasis is laid on 
disturbances of negligible curvature and practical methods of handling engineering 
problems arising in connection with the operation of locks or hydroelectric canals 
are given. 
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LIST OF SYMBOLS 

Symbols used in only one section are omitted unless they are of general signifi- 
cance. Important quantities in hydraulic equations are indicated as functions 
of x and t. All equations are valid in any consistent set of units. The significance 
of barred variables is explained in section II-l. 

A(x t t) cross-sectional area of channel. 

b width of free surface. 

exp ( ) exponential function. 

e base of natural logarithms. 

/( ), F( ) arbitrary functions, determined by boundarj^ conditions, 

/i substitution variable, (3co — 5Uo)gi/SUo(co — Uq). 

fi substitution variable, X /4 or — X /4. 

g acceleration of gravity. 

H(x,t) depth of water. 

Hq undisturbed depth of water. 

h(x,t) increase of depth due to disturbance. 

i slope of channel bed. 

k\ substitution variable, 1 — 3£/ /W 

ki substitution variable, 1 — Uq/ojq, 

I, m, n direction cosines. 

j) pressure. 

p a atmospheric pressure. 

Pxxi • • • generalized stress components. 

q discharge per unit width of channel. 

R xx , • • • turbulent stress components. 

R hydraulic radius. 

S xx , . . . viscous stress components. 

t time. 

U(x,t) mean velocity in channel cross section. 

Uq undisturbed mean velocity. 

u(x,t) increase in mean velocity due to disturbance (sections III— 3, ff.). 

u ^-component of velocity (mostly before section III— 3). 

u' velocity fluctuation due to turbulence. 

v ^-component of velocity. 

v' velocity fluctuation due to turbulence. 

w z-component of velocity. 

w' velocity fluctuation due to turbulence. 

x Cartesian coordinate, specifically, parallel to channel in bottom 

plane. 

y Cartesian coordinate, specifically, lateral to channel in a hori- 
zontal plane. 

z Cartesian coordinate, specifically, drawn upward with origin at 

bottom of channel. 

a Boussinesq coefficient of velocity distribution, 

e eddy viscosity. 

X coefficient of channel friction. 

\ Q coefficient of channel friction, H=H . 

fx viscosity. 

p density. 

a partial volume of wave per unit width. 

cr f partial volume of wave per unit width. 

r mean shear per unit area on channel wall. 

<P ( ) arbitrary function, 

x wetted perimeter. 

* ( )j ^ ( ) arbitrary functions. 

oj velocity of a wave-volume element. 

o?o velocity of a short wave, moving without change of form. 

I. INTRODUCTION 

This paper is the second of a series dealing with the motion of flood 
waves and other waves of translation in open channels. The purpose 
of the series is explained in the introduction to the first paper [l]. 1 

i Figures in brackets indicate the literature references at the end of this paper. 



Translation Waves in Open Channels 463 

Many of the significant properties of translation waves may be ob 
tained by assuming that the motion of fluid particles due to the wave 
disturbance is irrotational [1]. Experimental observation of these 
properties indicates that they are in good agreement with theoretical 
predictions, provided that the waves are moving in still water. If, 
on the other hand, the disturbances are produced in a liquid which 
is initially in motion, as in actual currents in an open channel, some 
disparity will be seen between observation and theory. Bazin's ex 
periments on the velocity of propagation of a solitary wave in still 
water or in a running current demonstrate these points [2], 

The divergence between observation and theory, when the wave 
is moving in a current, can be explained in two ways. First, irrota- 
tional or potential theory supposes that the original undisturbed 
motion is such that all the particles in a cross section move with the 
same constant velocity, U , which may of course be zero. This 
condition is not fulfilled for flow in actual channels, and it thus be- 
comes necessary to investigate the effect of the distribution of velocity 
in a cross section on the velocity of propagation of waves. An analysis 
of this effect has been carried through in great detail by Boussinesq [3]. 
With a view to simplifying the presentation and the solution of the 
problem, the authors give a new derivation (section III— 5). Al- 
though the results arrived at here are different in form from those of 
Boussinesq, it is believed that the difference has a qualitative rather 
than quantitative importance. 

This difference is due to two distinct forms of the equation of 
motion of unsteady flow in a channel, one adopted by Boussinesq and 
the other by the authors. The present method applies the theory of 
Reynolds stresses to the problem of the unsteady mean flow in a 
channel, without assuming any specific relation between the stresses 
and the velocity rates existing at points within the channel cross 
section. The effect of these stresses is evaluated in terms of the total 
friction at the wall, which is the important concept for mean flow in 
open channels and is usually available from observations on uniform 
flow. Boussinesq's method on the other hand is an extension of his 
own theory of turbulent flow, and specific relations between the 
turbulent stresses and the space-rates of the local velocities are as- 
sumed and introduced into the derivations. As the concept of 
Reynolds stresses is the starting point of the present analysis, the 
theory of these stresses is given in considerable detail (section II— 1). 
The equation of motion of a perfect fluid is the point of departure in 
interpreting the relation between the stresses and velocity fluctuations. 
This novel procedure simplifies the presentation considerably. It is 
hoped that the equation of motion for gradually varied unsteady flow 
in prismatic channels derived in this paper, eq 47, is sufficiently 
rigorous to merit the attention of river engineers. 

In the second place, the irrotational theory of translation waves 
ignores the effect of turbulence, which is one of the forces controlling 
the motion of the waves. Consider a short intumescence moving in a 
current. In the initial stages of motion the behavior of the wave is 
not affected by turbulence, since the wave is assumed to be short. 
When the wave has traveled a considerable distance the cumulative 
effect of friction becomes appreciable, particularly at the rear of the 
wave. In fact, experience shows that the front of the wave will 
travel with a velocity characteristic of irrotational waves, whereas at 
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the other end of the wave, the volume elements will move with a 
smaller velocity, the exact law of which is not as yet clearly understood. 
The wave thus flattens out considerably at the posterior end, with a 
consequent elongation. The tendency of short waves to deform in this 
manner is clearly demonstrated in Horton's experiments [4]. 

When the motion of the waves is controlled principally by friction, 
as in the quasipermanent regimes of rivers, the corresponding theory 
of translation waves is simple. But when it is supposed that the fric- 
tion and inertia effects are of like magnitude, analysis is confronted 
with a difficult task. The change of shape of a short intumescence 
during the entire journey of the wave in a current would be accurately 
described if a complete solution were available, which is not the case. 
However, some progress in this direction is possible if it is assumed that 
the effect of turbulence is small, requiring only a secondary correction. 
Such a theory has been worked out by Boussinesq, who assumed that 
the velocity of the wave-volume elements differs little from the theo- 
retical velocity of waves without friction [3] . The main purpose of this 
paper is to consider this theory and the corresponding solutions for 
waves of negligible curvature. In giving Boussinesq's solution, slight 
changes have been made, since it was desired to adopt Manning's law 
of friction in open channels [5]. Favre has employed the Boussinesq 
integrals to describe the motion of waves occurring during the opera- 
tion of locks [6] . Because of the engineering significance of the problem 
we have also reproduced these applications (section IV-9 and IV-12). 

The applications of the theory are presented in sections III and IV, 
which have been entitled first-order and second-order theories of 
wave propagation. The first-order theory neglects the vertical acceler- 
ation of the water particles and assumes that the combined effect of 
friction and velocity is small. It is found that this assumption requires 
that the waves be short. The second-order theory considers the effects 
of vertical acceleration but neglects the variations of the horizontal 
velocity components. As in Boussinesq's investigations, the concept 
of the velocity of a wave-volume element plays a very basic role. 
Equation 116 is the general expression for this velocity. Using this 
same formula, the authors have been able to study the gradual de- 
formation of a straight sloping front and also to compute the reduction 
of height of an abrupt wave front. 

II. FUNDAMENTAL PRINCIPLES 
1. REYNOLDS EQUATIONS OF MEAN MOTION 

In the Reynolds theory of turbulence [7; 8, p. 638; 9, p. 364; 10, 
p. 191] the velocity components u, v, w, in the directions x, y, z, 
respectively, are regarded as fluctuating functions of the variables 
x, y, 2, t. The fluctuations are frequent and rapid, and the periodicity 
of the_ changes is irregular. Denoting the mean values of u, v, w by 
u, v, w, and the mean value of the pressure, p> by p, we have for the 
instantaneous values 

u=u+u', 

v=v+v', 

w=w+w',} (1) 

and 

V=p+p', 
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where u', v', w' ', and f- are the fluctuations whose mean values u', v', 

w f , p f are equal to zero. The mean values may be temporal; that is, 
taken at a given point (x, y, z) over an interval of time of duration 
At; thus, for example, 



— i rt+M 

^ = -7-. I % 



udt. 



(2) 



Or again the mean values may be spatial; that is, taken at a given 
instant t over a space of volume AV; thus, for example 



l -AVJ Av 



udV. 



(3) 



For the present it will be supposed that the averages are temporal. 
A difficulty arises when taking temporal means in unsteady flow. In 
such cases it is necessary to make the time interval At long enough to 
include a great number of fluctuations, but short enough to exclude 

all but negligible variations of u, v, w, p. 

There are certain rules for taking the averages. Expressed sym- 
bolically [9, p. 366] 



and 



ab = 


^ab, 


a+b = 


=a+b, 


a= 


=a, 


da 


d — 



(4) 



where a and b are functions of x, y, 2, /; s represents any one of the 
independent variables, and a bar over a quantity represents an average 
of the kind just described. These rules are based upon assumptions 
that are not exactly applicable to the flow of liquids. The errors, 
however, are negligible when the fluctuations within the time interval 
At are sufficiently numerous. 

The specific contribution of Reynolds' analysis is the idea of in- 
terpreting the dynamic effects of the fluctuations u', v' , w\ as apparent 
tractional forces or stresses. The nature and the magnitude of these 
stresses can be readily obtained starting with the equations of motions 
of a perfect fluid. The dynamic equations in question are 



(5) 



du 

p Tt = 


-£+*' 


dv 

p ~d~r 


-%+*• 


dw 

p ~d~r 


-l+> z - 



where the operator d/dt, in terms of the local variations, is 

(15,5,5, 5 
dt dt dx dy dz 



(6) 
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Together with these, there is the equation of continuity, assuming 
that the liquid is incompressible, 



bu , bv jj ^_n 

da: by bz ~ 



(7) 



The liquid will always be assumed to be incompressible, in this paper. 
In eq 5, p is the pressure, p the density of liquid, and X, Y, Z the 
components of the force of gravity per unit mass of liquid in the 
direction of x, y, z. These relations are valid at any instant of time 
independent of the fact that the quantities, u, v, w, and p are fluc- 
tuating functions of time, provided that the velocities are not dis- 
continuous, and if discontinuities exist, they occur only occasionally 
during the time interval At. 

Replacing u, v, w in eq 7 by their respective values from eq 1, 
averaging each term, and usmg the last three relations in eq 4, together 
with the fact that u'=v'=w' = 0, we obtain 



bu _j_ bv , bw_ 
bx^by^bz"" ' 



(8) 



which states that the dilatation of the mean values of the velocity 
components vanish. It then follows from eq 1, 7, 8 that 



bu' .bv' _.bw' _ 
bx by bz ~ 



(9) 



The equations of motion, eq 5, are treated in a similar manner. 
We use 6, substitute for u, v, w their values from eq 1, multiply eq 9 
by pu' y pv' ', pw', respectively, and add to the left-hand side of each 
equation. Each term is then averaged in accordance with eq 4, 
assuming p to be constant, that is, that the fluid is homogeneous, 
obtaining 



du bp . v J)R X 
p dt == -&c +pX+ -te 



~ f "by~~ f ~bz~ 



bR X y ■ bRyy jbR Zy 



p dt~ by'* 1 * 1 T da; r by T bz 



dw_ bp. VJ bIt ztA bBy ZJ bR zz 
p dT-~bz +pZ+ W + by~ + ~te' 



(10) 



where, in this case, 



d b -b -b -b 

drbi +U b-x+ v dy +W bz' 



and 



R xx = — pu'u') R xy =R vx =—pu'v', 
Ryy=—pv'v'; R xg =R zx =—pu'w', 



R z , 



-pw'w' ' ; R vt —R zv = — pv f w' ., 



en) 



(12) 



We may refer to eq 10 as the Reynolds equations of motion for 
turbulent flow. 
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To interpret the above results we write the general equations of 
motion of a deformable continuum under the action of the external 
body forces X, Y, Z and surface tractions of an unspecified physical 
origin: 



dU — vxpPxxxpVvxxpVzx 

dv_ vipVxvypVw^pVzv 



r df 



dx by dz 



p ~dr pZ+ dx + by' + d^ ) 



(13) 



where the operator d/dt is the one defined in eq 6. 

It should be noted that these equations are merely Newton's law of 
motion as applied to any deformable continuum, and thus they apply 
to a wide range of phenomena. The left-hand members represent 
components of the rate of change of momentum per unit volume of the 
medium, whereas the right-hand members represent the components 
of the forces acting per unit volume. The nature of the medium 
determines the functional relation between the nine stress components, 

Pxx, Pvx, , and the other variables occurring in the equations 

and the intrinsic properties of the medium. Hydraulic equations 
established by integrating these equations over a cross section of the 
fluid normal to the flow are frequently referred to as momentum 
equations. 

The nine stress components completely determine the forces acting 
on a small surface dS at the point (z, y, z). If a normal is drawn to dS 
away from the body of which dS forms a part of the surface, and if I, 
m, n are the direction cosines of the normal to dS, then the x-, y-, z- 
components of the surface tractions on the body at the point (x, y, z) are 

Fx=lpxx+mpyz+np zx , ' 

F v =lPxv+mPw+np zy , (14) 

and 

Fz=lpzz+mpyz+np Z2 

Comparing eq 10 and 13, we see that 



Pxx=—P+Rx: 

Pyy=—P+Ry 
Pz.^—P+Rzi 



Pxy ~ Pyx = -Kxy ~ ^yx ) 
Pxz =z Pzz := H'xz =: tf'zx) 
Pzy = pyz —■ttzy = -ttyz ) 



(15) 



and thus the turbulent mean flow may be treated as ordinary stream- 



line flow where the stresses p XXj p T 



are identified with the 



stresses R xx , R xy These apparent stresses of turbulence are 

called Reynolds stresses, and represent mathematically the transport 
of momentum across a surface caused by the velocity fluctuations 
[10, p. 192]. 

If the flow instead of being turbulent is viscous in character, the 
stresses are of molecular origin. The stresses p XXi p zyi ... are then 
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identified with the viscous stresses S xx , S xy , ... in the following 
manner, 

Pxx :== - P\&zxj Pxy =:z Pyx =1 &xy z=: &yx) I 

Pyy=—P + Syy, PXZ=PZZ=SZZ = S ZX ; [ (16) 

Pzz=— jP+Su, Pyz=Pzy = S yz = S zy . J 

If the liquid is incompressible, these viscous stresses become 
du. Q Q /dvj)u\ 

o n ^ o o /dw .i>u\. /t ^ N 

where m is the viscosity of the liquid. 

It is perfectly feasible to carry through the derivation of the Eeyn- 
olds equations of motion, eq 10, using the viscous-flow interpretation 
of eq 13 as a starting point, and finally neglecting the viscous stresses 
in comparison with the Reynolds stresses, which is justifiable on the 
basis of experimental evidence [11, p. 186]. We have not adopted 
this method, since the results are identical with those obtained here 
by initially assuming a perfect liquid, that is, by neglecting the vis- 
cous stresses at the beginning rather than at the end, and since it is 
more complicated than the method used here. 

Later, in sections III— 1 and IV-1, we shall derive hydraulic equa- 
tions for mean flow in open channels, using eq 10 as a starting point. 
For this purpose, it is sufficient merely to suppose that R xx , R xy .... 
exist, there being no necessity for determining them as functions of the 
mean velocity gradients, because when the Reynolds equations of 
motion are integrated over the cross section of the channel, the 
quantities R xx , R xy , . . . either disappear or are expressed in terms 
of the wall friction, r , which is defined in connection with eq 45. 
Boussinesq, however, expressed R xx , R xy , ... in terms of the velocity 
gradients. He simulated eql7,_merely replacing S xx , S xy , ... by 
Rxx, Rxy, • • • ; u, v, w by u, v, w; and \x by e. The quantity e was 
called the coefficient of internal friction by Boussinesq [3, p. 46]; 
present usage seems to favor the term eddy viscosity, since e has the 
same dimensions as /x [11, p. 185]. 

When the channel is rectangular, Boussinesq assumed that e is 
constant for every point in a given cross section, but in a channel of 
constant roughness, directly proportional to the product of the wall 
velocity and the depth of the liquid. Such assumptions form the 
foundation for the developments in his "Essay." The validity of his 
theory is thus open to question, since the assumptions concerning e 
are widely at variance with present experimental knowledge. 

Having established the basic partial differential equations governing 
turbulent flow in open channels, we now turn our attention to the 
boundary conditions. Consider first the condition to be satisfied at a 
fixed solid wall, which we shall suppose is covered with uniform 
asperities of height k. A pillbox is erected on the surface, with height 
h, just sufficient to clear the asperities, and radius r, large with respect 
to 2h. (See fig. 1.) 
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Figure 1. — A pillbox enclosing the asperities on a channel wall. 

Since the liquid is incompressible, the net rate of flow across the 
surface of the pillbox must be zero. We neglect the discharge from 
the cylindrical side of the pillbox, since its area is small compared with 
that of the circular end. There must be at least one point on the end 
of the pillbox for which 

lu-\-mv J rnw=0 ) 

where /, m, n are the direction cosines of the outward-drawn normal. 
Averaging the terms in this equation and assuming that u, v, w vary 
slowly in an interval of length r, we have for all points on the fixed 
surface 



lu-{-mv-\-nw=0, 



(18) 



which is the desired boundar}^ condition for a fixed solid wall. If the 
wall is smooth, the argument is repeated, replacing the height of the 
asperities by the thickness of the laminar boundary layer. 

Next we consider a boundary condition which applies to the free 
surface. Let F (x,y,z,t)=0 be the equation of the free surface. As 
before in the case of the velocity components and the pressure (eq 1), 
the surface is considered to fluctuate rapidly and irregularly about 
a mean value, and we write 



F=F+F\ 



(19) 



where T is the mean value of F, and F' is the fluctuation whose tem- 
poral mean value F / is zero. 

Since particles on the surface remain there, we have the condition 

bF . bF , bF . bF n 
bt bx by bz ' 

where u s , v s , w s are the velocity components at the free surface. If 
we replace u 8 , v s , w 8 by expressions obtained from eq 1 and F by its 
value from eq 19 and average the terms, the result is 



bF. -bF . -bF .-bF A 



provided we neglect the terms u' s (dF'/bx) 9 



(20) 
These products 
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may be ignored, since either the correlation between u' SJ bF'/bx, 
etc., is zero or the term is small in comparison with those retained. 
On the basis of this assumption, the above equation is a boundary 
condition obtaining at the free surface of turbulent flow. 

In order to simplify the notation from this point on, u, v, w, p> 
etc., will be written for u, v, w } p, etc. Any further use of the bar to 
indicate an average will be in the macroscopic, rather than the 
microscopic, sense. 

2. EQUATION OF CONTINUITY FOR MEAN FLOW IN CHANNELS 

Up until this point we have derived equations which are valid for 
turbulent flow in general. We now derive an equation which applies 
to any tube of flow with a cross section which varies both in space and 
in time. This tube is assumed not to close on itself, and since it is a 
mathematical representation of a flowing channel, it is drawn as such 
in figure 2. The cross-sectional area of the channel, normal to x, will 
be represented by A=A(x,t), the average value of u over the cross 
section by U= U(x,t) . 




Figure 2. — Surface configurations at two different times in an open channel. 

The positive sense of x is chosen so as to make U positive. The 
actual direction of x is arbitrary, but for convenience it should approxi- 
mate the direction of the channel. 

The flow into the stretch Ax during the time At is 



ru+At 

i 



U{x^t)A{xr )} t)dt 1 
while the flow out of the same stretch is 



k 



U(x + Ax,t)A(x + Ax } t)dt. 
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Furthermore, the volume of liquid contained within the stretch at 
time t is 

A(x,t )dx, 

Jxo 

while the volume in the stretch at time t -\-At is 

rxo+Ax 

A(x,t + At)dx. 

Jxo 

Since the fluid is incompressible, net outflow plus increase of volume 
equals zero, or 

[U(xo+Ax,t)A(x +Ax } t)-U(x ,t)A(x ,t)]dt 

J to 



+ [A(x,t +At) -A(x,t )]dx=0. 

JXQ 



If we apply the theorem of the mean for integrals, divide by AxAt 
and then let Ax and At approach zero, we obtain 

b /TTA . \bA n (21) 

which is the equation of continuity for mean flow in channels. Since 
no assumptions have been made regarding the shape of the channel, 
this equation applies to open channels in all generality. 

III. FIRST-ORDER THEORY OF WAVE PROPAGATION 

1. THE EQUATION OF MEAN FLOW IN PRISMATIC CHANNELS 

Consider a prismatic channel and select as the origin of a right- 
handed system of Cartesian coordinates the lowest point of any con- 
venient cioss section (see fig. 3). Take the z-axis as the line parallel 
to the channel (not necessarily horizontal), with positive sense in the 



Figure 3. — System of coordinates in a prismatic channel. 

direction of flow, the y-axis as the horizontal line, and the 2-axis 
positive upward (not necessarily vertical). The channel slope, i, is 
the tangent of the angle 6 from the positive direction of the £-axis to 
the horizontal, measured in the counterclockwise direction, viewed in 
the direction of increasing y; in other words, a positive slope is down- 
hill. These conventions will be adhered to throughout the paper. 
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We assume that the angle of slope of the channel is so small that 
sin is sensibly equal to tan and cos to unity, hence X and Z in 
eq 10 are gi and — g, respectively. 

Again, in order to simplify the three equations of motion, eq 10, 
we make two distinct sets of assumptions, the first involving the 
terms on the left-hand side and the second involving the terms on the 
right-hand side. First, we assume that the second derivatives of the 
cross section, A, and of the mean velocity, U, and the squares and 
products of their first derivatives, may be neglected in comparison 
with the first derivatives. Such a flow is referred to by Boussinesq 
as gradually varied [3, p. 242]. In general, this requires that the 
sides of the channel be not much flatter than a 1:1 slope, that the slope 
of the wave profile be small compared to unity, and that the curvature 
of the wave profile be small compared to (1 / U) (bU / bx) or (1/U 2 ) 
(bU/bt). As a result of these restrictions, dv/dt and dw/dt are small 
with respect to dujdt, and the left-hand sides of the last two of eq 
10 vanish. 

Secondly, we make c ertain assumptions concerning the Reynolds 
stresses: (1) u / u / =v f v / =w / w f hence R xx =R yy =R 2Z ; (2) since u is the 
predominant velocity component, the derivative of the correlation 
between the fluctuations v', w' may be neglected in comparison with 
the derivatives of the correlations between w', v', and u', w', which is 
to say, that the derivatives of R zy may be neglected; (3) since the flow 
is gradually varied, the ^-derivatives of R xy and R xz may be neglected. 

On the basis of these two sets of assumptions, then, eq 10 reduce to 

bu . bu . bu . bu .lb, „ N . 1 d „ , 1 b „ ,„. 

bi +U bx + % + <%r^-p bx^- R - )+ -pby E ^ bz R "> (22 > 

0«^(p-B„), (23) 

and 

o=- p ^ 2 (p-R xx )+g. (24) 

The simultaneous solution of eq 23 and 24, with the dynamical 
boundary condition p—p a at z=H, is 

p—R zx =p a + P g(H—z), 

where H, a function of x and t only, is the ^-coordinate of the free 
surface, and p a is the atmospheric pressure. These symbols will be 
consistently used to refer to these quantities. According to the above, 
in gradually varied flow the pressure is hydrostatic everywhere. 
Differentiating the above equation, it is found that 

*^>-*£ M 

Each term in eq 22 is to be integrated over the cross section A at x 
In this way, the resulting equation of motion will involve only the 
mean flow, U, the cross-sectional area, A, and certain other quantities 
to be given later. We thus eliminate the variables y and z from con- 
sideration, since all of the quantities occurring in the final result will 
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be functions of x and t only. Such an equation is called a hydraulic 
equation, and bears the same relation to circuit equations in electricity 
as the general equations 10 do to the electromagnetic field equations. 
Let us first evaluate the integrals arising from the terms in the left- 
hand side. The mean flow, U, in a cross section, is implicitly defined 
by 



UA 



(26) 



= I udA. 

We write the above as a double integral 
UA= I dz I udy, 

JO J-Vx 

where y 2 and — jft are the ^/-coordinates of the two points on the wetted 
perimeter having the common distance z from the x, ?/-plane (see 
fig. 4). Obviously, y 2 and —y x are functions of z only, for a given 
prismatic channel configuration. Differentiating both sides of the 
above equation with respect to t, we obtain 



b (UA) 






where b 2 and —b\ are the ^/-coordinates of the two points on the wetted 
perimeter at the distance H from the x, ?/-plane. 




Figure 4. — Double integration with respect to y and z; cross section, looking upstream* 
524586—43 4 
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We now define u s to be the mean value of the surface velocities, 
u S) over the width b=b 2 -{-bi of the free surface: 

rbi 

u s b= u s dy, (27) 

J — b\ 

and the next to the last equation becomes 

since b(bH/bf) = 'dA/bt. Introducing the velocity coefficient a u implic- 
itly defined by 

ai U=^u 8 —U, (29) 

eq 28 may be written 

i>u, A A bU tt^A /r> _ N 

A i>t dA = A -dt- aiU w (30) 

We now turn our attention to the second term on the left of eq 22. 
We define another coefficient of velocity distribution a, as the ratio 
of the mean of the square of velocity to the square of its mean, both 
means being taken over the cross section, thus 



/ 



*U 2 A 



= (u 2 dA. (31) 



This coefficient is the only one which appears in the final equation, 
and was originally introduced by Boussinesq. It also plays an 
important role in the equation of steady flow in open channels [12] 
and must not be confused with the so-called Coriolis coefficient, which 
is frequently denoted by a. Assuming that a does not depend on x, 
and differentiating both sides with respect to x, in the same manner 
as was done in obtaining eq 28, we have 

Defining u a 2 to be the mean value of the squares of the surface veloc- 
ities over the width, 6, of the free surface, 

^?6= ( bl u*dy, (32) 

then the previous equation becomes 

since b(dH/bx)=(>A/dx. Introducing a third coefficient of velocity 
distribution a 2 , defined by 

' a 2 2U 2 =u s 2 -aU 2 , (33) 

we obtain from the last two relations, 

u s dA=aUA-^ aoU 2 -^-' (34) 

L 02 ox ox J 



h 



Next consider the integrals of the third and fourth terms of the 
left-hand side of eq 22. For the third term we write 



f du, A f* f* bu, 
J a 5y Jo J- yi by " 
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where the notation is the same as that used in the above two integrals 
(see fig. 4). The integral with respect to y is integrated by parts, 
obtaining 



JJ3H> 



Wo)2-M 



l]dz ~~L 



£>V 7 A 

a dy 



(35) 



where v Q and Uq are the velocity components at the wetted perimeter 
and the subscripts 1 and 2 indicate the two sides of the channel, that 
is, the points (—y h z) and (2/2,2), respectively. 




Figure 5. — Double integration with respect to z and y; cross section, looking upstream. 



Again, for the fourth term 



I w^-dA= ay \ 
JA <>Z J-61 J J Zi 



oz ' 



where Z\ is the ^-coordinate of a point on the wetted perimeter at the 
distance y from the z,2-plane (see fig. 5). This time the integral 
with respect to z is integrated by parts, obtaining 



w^rjAA 



bz* 



J-bx 



(w s u s —w u Q )dy- 



■I 






(36) 



where w , u are now evaluated at the point (?/,2i). 

If dx denotes an elementary arc of the wetted perimeter, x, measured 
in the counterclockwise direction, and the direction cosines of the 
inward-drawn normal are 0, m, n (figs. 4 and 5) then dy=ndx and 
dz——mdx ) andjeq 35 and 36 become 
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JM dA=± - f x ^md x - J A n^ y dA, (37) 

and 

w^-rdA— w s u s dy— J w u ndx— I wz—dA, (38) 

We add these two equations together, noting that as a result of eq 18 
the integrals with respect to x vanish, and combine the integrals in the 
extreme right by means of eq 8 and 34 obtaining 

Furthermore, since the equation of the free surface is 

F(x,y,z,t) = z-H(x,t)=0, 

and since bA/bt — b(bH/bt) and bA/bx=b(pII/bx) , the surface houndary 
condition, eq 20, becomes 

bA , bA 

Replacing w 8 in eq 39 by its value from the above equation, recalling 
that A is not a function of y, evaluating the integrals by means of 
eq 26 and 32 and eliminating bU/bx by means of eq 21, we have 

Introducing a final coefficient of velocity distribution, 

Ua^u-aU, (40) 

the previous equation becomes, using eq 33, 

) A V by dA+ j A w br^ U bt + ^ U bx- (41) 

We have now succeeded in evaluating all the terms on the left 
of eq 22 in terms of functions of x and t only, and proceed to the 
terms on the right. In view of eq 25, the first two terms become 
merely 

9t$/4^$'j!&^todA=M-9*&' (42) 

since H is not a function of y or z. 

The integral of the last two terms in eq 22 is transformed by Green's 
theorem into a line integral around the periphery of the cross section: 
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Expressing the line integral in terms of x as was done for eq 37 and 38, 

Comparison with eq 14 and 15 shows that the integrand of the line 
integral is the ^-component of the force which the fluid exerts on the 
channel walls and the free surface. If we denote this shearing force 
by r, and assume that the shear on the free surface is zero, and define 
r to be the mean shear of the fluid on the channel walls and bottom 



To 
then we have 



jjflfe, (43) 



Our choice of axes makes r positive, since r has the same direction as 
U. We now define X the channel coefficient of friction, 



j J2 
To^^—rT* (45) 



and eq 44 becomes 



m»« 



■fr-y ■ 



+Z-B„ Ua=-^- (46) 



Collecting the results of eq 30, 34, 41, 42, 46, and dividing them by 
Ay the desired integral of eq 22 is 






where R=A/x is the hydraulic radius of the channel. This is the 
desired hydraulic equation of motion for gradually varied unsteady 
flow in prismatic channels. It will be noted that only one coefficient 
of velocity distribution enters into the final equation, the one related 
to the distribution of the squares of velocities in a cross section. 

This equation was derived on the basis of the assumption that the 
flow was turbulent, which is the case usually met with in practice. 
Nevertheless, the equation also applies to viscous laminar flow, since 
the dynamical equations of motion, eq 13, are the same for both types 
of flow, except that in laminar flow the quantities j) XX} . . ., are re- 
placed by S X x, • • • ? rather than by R xx , .... It is tacitly assumed 
that the flow pattern is such that the assumptions made concerning 
the Reynolds stresses in deriving eq 22, 23, 24 are also true for the 
viscous stresses. Finally, in laminar flow a and X become quantities 
that may be calculated rather than ones that must be experimentally 
determined. 

2. COEFFICIENT OF RESISTANCE AS A FUNCTION OF H 

If it is assumed that the coefficient X defined by eq 45 is a function 
of H only, and is therefore independent of dU/dt and dU/dx, then X can 
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be evaluated from the ordinary parameters of uniform flow. Since, 
when the flow is uniform, V and H are constants, we have from eq 47 

. XUo 2 
or _ 

7&rV? (481 

where U Q and R , respectively, are the mean velocity and hydraulic 
radius, for uniform flow. 

We assume further that the propagation of the waves is restricted 
to that type of turbulent flow in which the friction is determined by 
the surface asperities independently of the Reynolds number. The 
resistance law is then represented adequately by Manning's formula 
in the following form [5, eq 61, 80]: 

where k s is the equivalent sand roughness [5, p. 713] of the channel 
surfaces. 

Equating the right-hand members of eq 48 and 49, we obtain 



\=0.0303 



/k\ 1/3 - 



Differentiating with respect to H 0) the value of H for undisturbed 
flow, we have 

d\ _ JK_dRo , ' 

dH<T 3i? dHo K0U) 

since 7? is a function of H . For a wide rectangular channel, R=H> 
and 

d\ _X^ f . 

dH " 3iJ # { } 

3. VELOCITY OF PROPAGATION OF SHORT WAVES WITH 
NEGLIGIBLE CURVATURE 

Assuming an initial condition of uniform flow in a very wide 
rectangular channel, let H be the depth of the liquid and U be the 
mean velocity of flow in a cross-section. If the flow is disturbed, this 
can be expressed by writing: 

H(x, t)=H +h(x, 0,1 
and (52) 

U(x, 0=l7o+u(x, *),J 

for the depth and for the mean velocity. The quantity u will, in 

general, henceforth signify the variation from Uq, the undisturbed 

value of U, and h the variation from H , the undisturbed value of H. 

We replace H and U in eq 47 by their values from eq 52. Since the 
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channel is rectangular, A is a linear function of H; since it is wide, 
B=H; and H and U are constants, hence eq 47 and 21 become 

bh .bu , , T7 , vdt*. /1 .Uo+udh . \(U +u) 2 /ro . 

^+dT +a(C7 °+^ +(1 - a) S^FI ST^2!CT' (53) 

We now introduce two restrictions of a different kind. First, we 
assume that a is very nearly equal to unity and hence the term con- 
taining I— a may be neglected. Secondly, we assume that 

(/) ^<^, (55) 

en) rf-tir*??' ( 56 ) 

(///) , "5i + *5 < 5r (57) 

On this basis eq 53 and 54 reduce to 

We differentiate eq 58 with respect to # and elinunate the terms 
containing u by means of eq 59. This gives 

A particular solution of this is 

h=F(x-a> t), (61) 

where F is an arbitrary function, provided w is a root of 

u 2 -2U o>o+U 2 -gH =0. 
The two roots are 

co =^±V<75o; (62) 

and it is evident upon inspecting eq 61 that the wave is propagated 
without change of form with the celerity, or velocity of propagation, 
co in the positive z-direction, since at time t-\- (Az/coq), the wave profile 
is displaced by Ax from its position at time t. When an observer is 
moving with the current, that is, moving with the velocity U Q) a wave 
may be moving either downstream or upstream with respect to the 
observer. In the first case, C/ <C w o, and the wave is descending; in 
the second case, £7o>«o> and the wave is ascending. Accordingly, 
the larger root of u Q gives the celerity of a descending wave and the 
smaller root the celerity of an ascending wave. In all formulas 
containing an ambiguous sign the upper sign will pertain to descending 
waves and the lower sign to ascending waves. 
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It is now a simple matter to determine the function u. From eq 
61, we have 



bh_ bh 



(63) 



Substituting this value of bh/bt in eq 59, it becomes 
whose integral is 



bu Uq—o) Q bh___~ 
bi + ~~H~ bx~ V ' 



u+U^h+f(t) = 0. 

Since for any time t, at A=0, ^=0, then since the liquid is undisturbed 
at infinity /(£)=0, and 

u=(u -Uo)^, (64) 

or in view of eq 62 

which states that the excess of the mean velocity of the particles in a 
cross section under the wave is proportional to the wave height A. 

It will be noted that the above equation for the celerity of waves, 
eq 62, is the same as that developed for the motion of irrotational 
waves in still water [1, eq 56], except that a constant velocity of flow, 
Do, is superposed on the wave velocities relative to the water, -y/gH . 
This would appear to indicate that the friction due to turbulence 
which is inherent in the undisturbed flow, t/ , has no effect on the 
propagation of waves. In general, the effect of turbulence cannot be 
neglected, but in this case the assumptions made in deriving eq 62 
restrict the size of the wave in such a way that the effect of turbulence 
is negligible. 

The physical significance of these mathematical assumptions will 
now be investigated. Consider first 

/ttt\ bh , ibu ^bh 

(III) %r+A5r-<5^" 
bx bx bt 

Replacing u hy its value from eq 64, and bh/bt by its value from eq 63 ' 
we have 

h 



H<<-2\ 



Consider secondly 

/t\ bu bu 

(I) V^aT 

Using the same relations as before, we have 



H Uq—C0 

Obviously, the first condition is sufficient for the second. 
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Botli of these restrictions have no bearing on the fact that the 
undisturbed flow is controlled by friction. They are necessary but 
are not sufficient conditions for a wave to move without change of 
form. A more significant restriction arises from the assumption that 
the right-hand member of eq 53 is negligible: 

(IL gi "2(H +h) <<: ar 

Expressing the left-hand side in terms of u[U and h/Ho, and replacing 
X by its value from eq 48, we obtain, since JS =H , 



i/, o uH /uVII \ 



Eliminating u via eq 64, and transforming the right-hand side by 
means of eq 63, we have 



^■0- %-fe-O'a) ^ 



-co (co — U )l 1+Jj) 



h \ \H J dx 



In contrast to conditions (I) and (II), where we obtained an 
inequality comparing A/Ho with a constant term which is determined 
by the characteristics of the channel, we find in this case that both 
sides of the inequality are functions of h. We are also confronted with 
the fact that bh/dx generally vanishes at the apex of the wave profile, 
and hence the inequality cannot be satisfied at that point even if the 
left-hand member also vanishes. 

In order to arrive at a workable criterion for determining whether 
the theory is applicable, we integrate both sides of the inequality 
with respect to x. thus obtaining an average condition: 



jgHo 



-co (co — U ) 






bh\ 

dx\ 



dx. (67) 



It is not feasible to integrate the left-hand side of this inequality unless 
the function h(x,t) is given. It is therefore necessary to assume that 



™ (E- 1 )!.* 3 -!?' (68 > 



and 



and eq 67 reduces to 



(V) £r<l, (69) 



4feKfSh (7o) 
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where 

co (7o(co — (Jo) 

If the wave profile is such that (IV) and (V) are not both satisfied, 
then eq 70 does not apply, and the integral on the left of eq 67 must 
be evaluated. Assuming that the wave profile has a unique maxi- 
mum height h m , that bh/dx is everywhere finite, that L is the length 
of the wave, and that h is always positive, then eq 70 becomes 

L 2 

H «W (71) 

where j is a fraction such that the volume of the wave per unit width 
of channel above the undisturbed surface is jh m L. This last inequality 
is the desired average condition, one which restricts the length of the 
wave. Waves fulfilling this condition will be referred to as short 
waves. 

In addition to this restriction on the length of the wave, we still 
have the three conditions on the maximum height of the wave, 
expressed by eq 66, 68 and 69; the first being derived from the ori- 
ginal restriction (III), and the other two from the approximations 
used in obtaining eq 71. Each of the first two restrictions states 
that the maximum height of the wave must be small with respect to 
a constant determined by the channel parameters. For convenience 
we shall express these constants in terms of U /-yJgH . 

For descending waves, u =U J r'y/gH Q , and eq 66, 68, 69 become 



and 



(IIIa) So«-KJro +1 > 
h * u ° ( & A 



(IVa) ^« Uo 



(V) ^-«1. 



For ascending waves, co =£7 — -y/gH , 
and the same three equations become 



^(vSrO 



a 



and 



(IVb > A«S(vtr„+ 2 ) 



(V) j|«l. 



These relations are made more intelligible by means of figures 6 
and 7. The quantity h/H Q must be small with respect to all of the 
quantities plotted in the figure, that is, small with respect to the 
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least of them. It will be noted that the theory breaks down com- 
pletely for ascending waves in water of critical depth, that is, it does 
not apply to the case of the hydraulic jump. The fact that it breaks 




Figure 6. — Limitations on height of short descending waves moving without change 

of form. 
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Figure 7. — Limitations on height of short ascending waves moving without change 

of form. 

down for Uo/^gHo=2 wiien the wave is descending is due merely to 
the approximation used in deriving eq 71, and is not inherent in the 
theory. 
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Our three original assumptions thus indicate that if a wave is short 
and of small height, and moving in a wide rectangular channel where 
the undisturbed current is controlled by friction, the waves move 
without change of form, and their velocity is constant. It was also 
assumed in deriving eq 47 that the surface curvature is negligible, 
and this restriction is still in force. The simple theory thus gradually 
loses its applicability as the waves become larger and larger. This 
case will be treated in the later sections. 



4. PROPAGATION OF DISCONTINUITIES OF SLOPE 

One of the simplest types of wave profiles which we can imagine 
superposed on the initial uniform flow in a wide rectangular channel 
is one which is polygonal in form, such as a triangle. We shall inves- 
tigate the motion at a vertex of such a polygon, where the slope of 
the wave is discontinuous. 

We assume that at t=t Q , the discontinuity of the slope bh/bx is 
located at the point x—x , and the water surface is represented by 
two straight lines (see fig. 8). The mean velocity, U, and the depth, 
H, will each be described by two distinct functions; downstream from 
x , we write U=Ui and H—Hi, and upstream from x , we write 
U= U 2 and H=H 2 . We assume that both Hi and H 2 have negligible 
curvature, and that the channel is rectangular. Under these condi- 
tions, eq 47 applies to each pair of functions, and we write, putting 
a equal to unity: 



bx 



bt 



r^+^r+Ui 



bH 2 bU 2 

dx + bt 



g 



■u 2 



bUi 

bx " 

bU, 

bx 



. XJ7, 2 



--gi- 



2H 2 



(72) 



(73) 




+ 4f 



Figure 8. — Motion of a discontinuity in surface slope. 

In a similar manner, we obtain from eq 21, since A is proportional 
tofl, 






dt 



dx 



dx' 



(74 
(75) 
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Following Hadamard [13, p. 83], we introduce the following nota- 
tion for the change in value of a derivative as the discontinuity is 
traversed : 



i = R>#1 dff 2 dt/j TdET] 



ag 2 dg!_rafli bU 2 i>U l _ 

dx dx 



Subtracting eq 73 from eq 72 and eq 75 from eq 74, we have, since 
H, U, and X are continuous, 

where H and U are evaluated at x— x Q , t=t . Since four disconti- 
nuities are involved, we need two more relations. 

In figure 8 the displacement of the discontinuity during the time 
interval At is indicated. If AU\ and AU 2 are the changes in the 
mean velocities, then 

ATT * U U . ^a, ( 78 ) 

ATT * U *A J_ & ^A, (79) 

Similarly, if AHi and AH 2 are the changes in the height, then 

«-£«•+$» <81) 

Subtracting eq 79 from 78 and 81 from 80, 

KCM3£> 
°=[f>+[f>. 

since C7 and H are continuous. We denote the velocity of propagation 
of the discontinuity, Axj\t> by 0. Thus the last two equations become 
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We now have a system of four homogeneous linear equations, eq 
76, 77, 82, 83, in the four unknown discontinuities. Since the dis- 
continuities do not vanish, the determinant of the coefficients must 
vanish, and this relation gives 



Q-U=±J g H. 



(84) 



The velocity at the front of a descending wave, as shown in figure 
9, will now be considered. In this case, 2 is positive, and we see from 
eq 83 that [bH/bx] and [bH/bt] have opposite signs, and since in front 
of the wave, bH 1 /bt=0 ) bH 1 /bx=0, then behind the front of the wave, 
bH 2 /bt^>0, bH 2 /bx<^0. Now, U= U and H=H at the discontinuity, 
and the velocity of the wave^front is 



0=Z7 +VfffloJ 



(85) 



the positive sign is selected, since O>0, regardless of the magnitude 
of U and H Q . We thus see from eq 62 that at the front of such a 
wave Q = o) . This analysis fails to give any information concerning 
the change in the magnitude of the discontinuity. Such information 
may be obtained by using Boussinesq's equation for the velocity of 
propagation of a wave-volume element. A typical case is discussed 
in section IV-6, page 498. 




Figuee 9. — The front of a descending wave, 

5. EFFECT OF VELOCITY DISTRIBUTION ON THE CELERITY OF 

SHORT WAVES 

In the derivation of the celerity of short waves, given in section 
III-3, page 478, it was assumed that 1 — a is a small fraction that 
may be neglected. This is generally true in practice, but we now 
drop this restriction and thereby consider the effect of the velocity 
distribution on the celerity of short waves. Retaining a and using 
the assumptions of section III-3, eq 55, 56, 69, the equation of motion, 
eq 53 now becomes 



bh bu , TJ bu ( 

g b-x + bi +aU «-bx +(l ' 



N U +u bh A 



(86) 



Even though 1 — a: is not negligible, it is still true that in practice 
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a is of the same order of magnitude as unity. From eq 69 we 
have therefore 

n_ x cop— U dh h ^p— Up dh 

u a) II dt W<T Ho M' 

Using the relation between h and u obtained in the previous analysis, 
eq 64, this inequality reduces to 



„ , u dh du 



and eq 86 thus becomes 



Differentiating with respect to x, assuming that a does not vary 
with x, and eliminating the terms containing u by means of eq 59, we 
find 



b 2 h . r T ~b 2 h , f TT2 zr^(> 2 h n 

W+ 2aUo mt +(aUo '-0®$?-*' 

A particular solution of this is, as before, 

h=F(x—u t), 
where F is an arbitrary function, provided co is one of the roots of 

a> Q 2 -2aU o>o+(aU Q 2 -gH ) = 0. 
The two roots are 

^o-(H-r7)^o±V^o+(^ + ^ 2 )C/o 2 , (87) 

where 97 =a— 1. (88) 

This is the formula giving the celerity of short waves when the 
effect of the velocity distribution is considered. It will be seen that 
rj is the excess over unity of the ratio of the mean of the velocities 
squared, taken over the different points of a section, to the square of 
the mean velocity. A corresponding expression is given by Bous- 
sinesq [3, p. 285, eq 265], which differs from the above only in the 
values of the coefficients. 

It is preferable that the value of 7? for a given case be ascertained by 
observation. If this is not feasible, a value may be obtained by 
considering the relation which exists between rj and U Q and u m , the 
latter being the maximum velocity occurring at the surface of a wide 
channel. A relation of this kind can be derived if it is assumed that 
in a channel of great width the "velocity defect/' that is, the dif- 
ference between the maximum velocity u m and the velocity u at a point 
whose height above the bottom is 2, depends only on the so-called 
"shear velocity" [5, p. 709] and the relative depth of the water, and 
is independent of the character of the asperities on the channel 
bottom. Mathematically speaking, 

u m —u=u*j{{), (89) 

where ^ = z/H 0) and u* is the shear velocity, 

v*=JwB*. (90) 
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P Multiplying both sides of eq 89 by d$ } integrating from z=0 to 
z=H that is, from to 1, putting 

and recalling that the mean value of u is U 0) we obtain 

u m -U Q =NiU*. (91) 

Squaring both sides of eq 89, we obtain in a similar manner, using 
eq 88, 31, 26, 

uJ-2u m U +(l+ v )U 2 =N 2 u* 2 , (92) 

where 

J j({) 2 d£=N 2 . Eliminating the shear velocity, u*> between eq 91 
and 92, and solving for 77, we obtain 

which is the desired relation expressing 77 in terms of u m and H7 . 

To determine the numerical values of N\ and N 2 we make use of 
some of Nikuradse's experimental results concerning the velocity dis- 
tribution in sand-coated circular pipes [14, p. 18-19], The data in 
table 1 were obtained by selecting three determinations of the velocity 
distribution for each of six relative roughnesses. In Nikuradse's 
paper the data for each relative roughness are given in a separate 
table, and the three determinations on the extreme right of each table 
were selected, except in Zahlentafel (table) 13, where the first one and 
the last two were chosen. In each case, velocity-defect ratios were 
computed for each value of y/r, which corresponds to z/H in our 
notation, and the mean of these 18 determinations was entered in 
table 1. Plotting curves for /(f) and [/(f)] 2 from these data and in- 
tegrating numerically, it is found that in circular pipes, iVi=2.52 and 
N 2 = 12.29. Since the "velocity defect 7 ' relation is the same both for 
circular pipes and for open channels of great width these values of 
Ni and N 2 may be substituted in eq 93, which becomes 

,=0.94(^-1)*. (94) 

An expression for r? in terms of X may also be derived. Eliminating 
u m between eq 91 and 92, we have 

Eliminating u* and U via eq 90 and 48: 
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which becomes, when we replace N 2 and N\ by their numerical values 

77 = 2.97\. (95) 

Table 1. — Velocity defect as a function of distance from channel bottom. 
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It will be recalled that in deriving eq 62 of section 1-3 it was 
assumed that a is nearly unity, that is, that n is small, and that this 
restriction was removed in deriving eq 87. By comparing these two 
equations we can establish an approximate relation;, which will show 
how small rj must be in order that eq 62 may be used. j g 

Equation 87 may be approximated by 

co =(l + ^)C/o±V^o(l+^- 2 JD 
and this reduces to eq 62 if 



iU ± 



v+v' 



JgH 



<U ±^/gH Q . 



Dividing by U , neglecting r? with respect to rj, using eq 48, this 
becomes 



«wm 



a 



or, replacing i? by its value from eq 95, 



2.97X 



Nvfh+V 



2% 



(96) 



(97) 



Bazin's inexhaustible store of experimental data furnishes an 
example [15] of the improvement obtained by using eq 87 instead 
of eq 62. The observed values, together with the results of the two 
formulas, eq 62 and eq 87 are shown in table 2. These experiments 
were conducted in the discharge channel of the Grosbois reservoir. 
This channel was constructed of masonry, with the steep slopes of 
0.101 in the upper part and 0.037 in the lower part. Before entering 
the channel, the water drawn from the reservoir flowed over an apron 
11.20 m wide just below the gates, and which narrowed down to 1.80 
m (the width of the channel). The axis of the apron described an 
arc of 90°, with a radius of 40.50 m, in order to join with the channel. 
As a result of the narrowing of the cross section and the change of direc- 
tion of the apron, the flowing water was subjected to various reflec- 
tions at the walls and these complex movements set up waves which 
passed rapidly down the channel at fairly regular intervals of about 2 

524586—43 5 
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sec (31 per minute). These waves were rather small in height, but 
their transit was easily observed owing to the slight noise which 
accompanied them. 

\ The observed values of w given in the table 2 are the averages of 
six individual observations, the ± term representing the standard 
deviation. The values of rj are computed from eq 94. The table 
shows that the values of co computed by the simple formula, eq 62, 
are considerably in error, whereas the values of co obtained from 
eq 87 compare reasonably well with the experimental values. 

JH Table 2. — Experimental and calculated values of coo for torrential flow 



Experimental data (^=9.81 m/sec 2 ) 


i 


Ho 


l/o 


Mm 


COO 


0.101 
.037 
.037 


m 

0. 110 
.150 
.235 


m/sec 
3.785 
2.744 
3.481 


m/sec 
5.51 
3.49 
4.55 


m/sec 
6. 25 ±.07 
4. 32 ±.08 
5. 75 ±. 06 


Computed values eq 62 


Computed values eq 87 


coo 


Error 


V 


coo 


Error 


m/sec 
4.82 
3.96 
5.00 


% 
-22.9 
-8.3 
-13.0 


0.195 
.070 
.089 


m/sec 
6.63 
4.36 
5.65 


% 

+6.1 
+0.9 
-1.7 



IV. SECOND-ORDER THEORY OF WAVE PROPAGATION 



1. EQUATION OF MEAN FLOW WITH APPRECIABLE VERTICAL 

ACCELERATION 

Assuming that the liquid is flowing in a wide rectangular channel 
(t?=0, and R VX =Q)) that the variation with z of the ^-component, u, 
of the local velocity is negligible, or u= U(x, t) and that the assump- 
tions made in section III— 1 regarding the Reynolds stresses and the 
orientation of the axes still hold, then the Reynolds equations for 
turbulent flow, eq 10 become 



(>t^ U dx~ 



--gi- 



1 b(p-R xx ) 
p dx 

b(p—R xx ) 



1 dR z 



dz 



(98) 



and 



by 



bw . j-jbw . bw 1 b(p—R xx ) 

bt + U bx +w bz~ ~ P bz 



-9- 



(99) 



We shall eliminate p from the eq 98 with the aid of eq 99, and then 
eliminate z by integrating over the cross section. 
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We begin by evaluating w in terms of the other variables. Since 
u= U and v=0, the condition of continuity eq 8 reduces to 

bw = dff 

dz dx 

Integrating with respect to z % we have, since w=0 when 2=0, 

dU 

w= — z^> 
ox 

which, using a condition of continuity of mean flow, derived from eq 21 
on the basis of the assumption that A is proportional to H, 

reduces to 



w~- 



=KW+«f) tm 



Replacing w in eq 99 by its value from eq 101 and neglecting squares 
and products of the first derivatives of H and U compared with the 
second derivative, we obtain, since H and U are functions of x and 
t only, 

d(p-R xx ) _ P zU 2 /d 2 H 2 d 2 H 1 d 2 H \ 

dz ~ pg H \bx 2 + Udxdt + U 2 dt 2 / 

Multiplying by dz and integrating, 

U 2 HS z 2 \fd 2 H , 2 yff , 1 d 2 



V-Rxz=Pa + Pg(tt-z)+P^{l-Jj 2 )y-fa2 + 



UdxdVU 2 dt 



since at z=H, p—R xx =p a ' We now differentiate the above equation 
with respect to x, again neglecting the products of the lower deriva- 
tives as compared with the higher derivatives, and the value thus 
obtained for d(p—R rx )/dx is inserted in eq 98. The resulting equation 
is multiplied by dz and integrated from z=0 to z=H. Since (R xz )z=o 
equals r , the shear on the channel bottom, and (R xz ) z =h is zero, we get 

<>E, T7 dU dff , HIP/&H ,2 d»g . 1 d 3 H\_ . \U 2 , . 
dJ+^dz^dz^ 3 V &? r+ U dx 2 dt ^ U 2 dxdf)~ g% 2H J {WZ) 

upon replacing r by its value from eq 45. 

This is the equation of mean flow for a wide rectangular channel 
where the vertical accelerations are not negligible but the horizontal 
velocities in a section are practically uniform. 

2. HEIGHT OF WAVE AS A FUNCTION OF x AND t 

Replacing dU/dx in eq 102 by its value obtained by solving eq 100, 
and multiplying through by H, we obtain 

, 2 yg , 1 b 3 H\ „. x rr2 /in _, 
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We again introduce the expressions 

H(x, t)~H +h(x, t),\ 
and [ (52) 

U(x, t) = U +u(x, t)\ 

Substituting in eq 100, we have 

bh . Tj-bu TT bh ,c)u dh 
M +H °bx+ U °bx +h dx +U ¥x =0 ' < 104 ) 

while eq 103 becomes 

TT bu TT bh , t TT TT2\^h , idu dh ■ . / , orT N bA 

+-3~W + )K235i + I7? bxd?) =9Hl ~2 U > (105) 

when ^ 2 is neglected with respect to t7 2 , and in the term containing the 
third derivatives, w and h are neglected in comparison with V and £T. 
Various assumptions as to the nature of the right-hand member will 
be made in later sections. 

The terms of small order in the left-hand members of the above 
equations, that is, the last two in eq 104 and the last four in eq 105, 
are simplified by introducing values from the first approximate solu- 
tion: 

and 

«—*&• < 63) 

From the latter and eq 61, we obtain the further relation 

m_ 2 d% 

For the sake of brevity, we introduce the new variables 

3EV 



and (106) 

o> 
and we also make use of the equation 

gH Q ={^-U,)\ (107) 

which is derived from eq 62. On the basis of the above relations, then 
eq 104 and 105 become 
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5i+<+<+<--<(I)=°; W 

and 

We now differentiate eq 108 with respect to t. In differentiating 
the last term, which is of the second order, we use the relation 

since these small terms may be treated as functions of x—o) t [1, p. 73]. 
We next differentiate eq 109 with respect to x, and subtract this equa- 
tion from the one obtained from eq 108. The result is 

which is the fundamental equation for determining h as a function of 
x and t. 

This equation, which was derived from eq 102, has the same range 
of validity, namely, it applies to flow in wide rectangular channels 
where the horizontal velocities in a section are practically uniform. 
This equation takes account of vertical accelerations, finite wave 
heights, surface curvature, and friction. 

3. CELERITY OF A WAVE- VOLUME ELEMENT 

In order to integrate eq 111, the concept of celerity oj a ivave-volume 
element, an essential simplification due to Boussinesq [3, p. 451], is 
introduced. To define this concept, we consider the volume of the 
wave per unit width of channel, above the undisturbed primitive level 
of the liquid and between x, an arbitrary point, and oo f for a descending 
wave, and between x and — °°, for an ascending wave. Denoting 
this volume by o-, we have 



= I hdx and a= hdx, 



for descending and ascending waves, respectively. Recalling the 
convention of signs introduced in section III— 3, the above equations 
become 

'"hdx. (112) 



-c 



Imagine a plane, normal to the channel axis at the point x f and which 
moves in such a way that the wave volume in front of it, a, remains 
constant. This velocity is the celerity of the volume element <r, and 
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is denoted by a>=oo(Xj t). Differentiating eq 112 with respect to t, 
we have 



r ±co bh 

0=± Jx bt dxTha > 



since by definition dx/dt=w, a remains constant, and we assume that 
the surface of the liquid is undisturbed at infinity, that is, h vanishes 
at ±00. Differentiating the above equation with respect to x y we 
have, regardless of choice of sign, 

»+!*•>-* (H3) 

which gives the relation obtaining between w and h. 

We now introduce the value of dh/dt from eq 113 into the first two 
terms of eq 111, thus obtaining 

- Uo ( aa -u )^- Y -^+- T - w) = -dxK gm -2 u / 

This equation is immediately integrable once with respect to x. The 
arbitrary function of integration vanishes as a result of the previous 
assumption that the surface of the liquid is undisturbed at infinity. 
We thus have 

|(M +2 u£(M + (gH - W)~ 

We next transform eq 114 into a form which can be integrated. 
In this process it is necessary to operate only on the first three terms 
in the left-hand member. If we assume that co/w is near to unity, 
we may write as an approximation 

w =co (l+5(x,0). (115) 

Substituting this value of co in the terms under consideration, we obtain 

Now using eq 113 and 107 we obtain 

— ^o^(M + co ^(/^) + 2o;o?7o-^(^ + ^) + (coo 2 — 2co C7 )-^.' 

Applying eq 115 again to the first term and eq 110 to the second term 
(which can be done because hd is a small quantity), and cancelling 
terms, we obtain 

-2u ^ x (h8) + 2a> U ^(h6), 
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which can be written 

— u Q (u —U )^(2h5)' 

Hence eq 114 transforms into 
— o> («o — U ) gj (2*5) + co (co — J7 ) gj 

which, when integrated, gives 

4 H + 6ft dx 2 2A«o(«o-^o)J±.V 2^ /'*• 



and hence from eq 115 

a) . , 2 + fr, ft , hH 2 b 2 h , 1 f*Y „. \ TT2 \j , 11A . 



This is a general expression for determining the effects of wave 
height, curvature of the wave profile, and friction, on the celerity of a 
wave-volume element. The integral containing the friction term 
requires further consideration before the actual wave profile can be 
determined as a function of time. 

It remains to evaluate that part u of the velocity U under the wave 
in eq 52. Transforming the equation of continuity, eq 100, by means 
of eq 52, 

^ t +-^ x (U +u)(H +h)=0. 

The result of subtracting eq 113 from the above is 

^[(U +u)(H +h)-ho>]=0. 

Integrating, and using the boundary condition at x= co f for all times 
t } h=0, and u=Q, we have 

(U +u)(H +h)-hu>=U H , 
or 

-ggL (117) 

4. DEFORMATION OF THE WAVE PROFILE, ASSUMING 

A natural assumption, in proceeding from eq 116, is that the right- 
hand member (gHi—\U 2 /2) is negligible. From eq 56 this implies that 
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the waves are short, as described in section III-3. On this basis, eq 
116 becomes 

W= T + "TH + 1TW (118) 

or when we introduce the values of k u k 2 , and w from eq 106 and 62, 

— Di-±VOT(l+S+Sf-g> (119) 

Comparison with an earlier paper [1, eq 112] shows that this expres- 
sion is equivalent to that for the propagation of an irrotational wave 
relative to the mean velocity of flow, U . The assumption that 
gHi—\U 2 /2 is equal to zero is thus tantamount to assuming that the 
only effect of friction and slope is to superpose a constant velocity U 
on the wave, which itself is propagated without any change of shape 
due to frictional damping or to the x-component of gravity. 

The deformation of the wave profile, in a short time interval, is 
determined by using eq 119. Equation 113 can be written 

bh bh. d,, N 

For small values of £, the integral of the above equation is 

h=f( x -a t)-t^ x [h(a>-u Q )]. (120) 

In verifying this solution, we must use eq 110 when differentiating the 
second term on the right, which necessitates assuming that this term 
is small compared with the first term. In other words, the magni- 
tude of the time interval in which this approximation is permissible 
varies inversely with d[A(co— <a )]/bx. Clearly, since at £=0, we have 
h=f(x), the function j represents the wave profile at the fiducial 
instant. Finally, using eq 119 and 62, eq 120 can be thrown into the 
form 

Wfr-rf)TiVi^g+^g> (121) 

For a discussion of the deformation of a wave of negligible curvature, 
propagated without damping, the reader is referred to the earlier paper 
[1, p. 89]. All the formulas on the page cited are directl y ap plicable, 
provided the expression -y/gH(l + . . .) is replaced by U ±-JgH (l-\-. . .), 
wherever it occurs. This, again, is merely equivalent to superpos- 
ing the velocity of flow U 0j on the velocity of propagation of the 
wave-volume element. 

5. EFFECT OF FRICTION AND SLOPE OF BED ON THE PROPAGATION 

OF WAVES 

In the above solution of eq 116, the term gHi—\U 2 /2 was ignored. 
If the wave is short, the influence of this term is negligible, as explained 
in section III— 3. For long waves, however, these effects become 
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appreciable at distances sufficiently removed from the head of the 
wave. The problem of determining these effects is considered sepa- 
rately for the two cases, (I) ?>0 and C/ >0, and (II) i=0 and U =0. 
The first case will now be considered. 

From eq 52, we have, after neglecting u 2 in comparison with U 2 , 

gHi-~U 2 ^gHd-lu 2 +ghi-\U u. (122) 

Considering the first two terms of the Taylor expansion of \=\(H) 
at H Q) / /7\ \ 

x=x °+(Sr' (123) 

where \ =\(H Q )j and dX/dH is evaluated at FI=H . Replacing X by 
the above value, and neglecting the small hu term, eq 122 becomes 

gHi-\w=gI-I i-^UJ+ghi-\^hW-\TJ u. 

From eq 48 we have 

gH i-h U( 2 ==z0j (124) 

that is, in uniform flow the effects of slope and friction are in equilib- 
rium, and thus 

gHi- -U 2 =ghi~2 jjjhU 2 —\>U u. 

Evaluating d\/dH by means of eq 51 and 124, and introducing the 
value of u from the first approximation, eq 64, the above equation 
becomes 

gHi-\u 2 =-2{^-U )hj u (125) 

where 

Since from eq 112, da/dx=Th, eq 125 can be written 

gH i-\u>=±2^ a -U )l(^} (127) 

Substituting this in eq 116, we have 

or introducing the values of k x and k 2 from eq 106, and using eq 62 

„- Pl=± (VSH( 1+ ^ + §'g)-£) (129) 
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This integral of eq 116 is due to Boussinesq [3, p. 451], except that the 
value of /i given above is new, being based on Manning's quadratic 
law of resistance developed in section III-2. 

We now consider the sign of the quantity afi/h. If the intumescence 
is everywhere positive or everywhere negative; that is, for all x, h^ 
or h ^ 0, then a has the same sign as h, and tr/h is positive. At the head 
of the wave this ratio vanishes, since a is of higher order than h, and it 
increases toward the rear of the wave. Since/i is constant for a given 
channel and depth of water in uniform flow, the absolute value of the 
term varies in the manner just described, and its sign is the same as that 
of /i. Upon examining eq 126 we find that /j is positive if (1) £7o]>ca& 
or (2) Uq<^03 and 3co >5£7o, and negative if (3) Z7 <C^o and 3co <5?7o- 
Comparison with eq 62 shows that (1) applies to all descending waves, 
(2) to descending waves where U 2 <^(^/i)gH , (3) to descending waves 
where t7 2 >(9/4)#i7 . Thus, except in the last case, the combined 
effect of friction and slope on the damping of the wave is to cause a 
diminution of the absolute value of the velocity of propagation with 
respect to the velocity of flow of the various parts of the intumescence, 
the effect increasing with increasing distance from the head of the 
wave. The last case is of special interest. It has been shown else- 
where [16] that if [/ 2 >(9/4)^ff , the waves become unstable, a 
phenomenon which is apparently closely related to the origin of roll 
waves. 

6. DEFORMATION OF A STRAIGHT SLOPING FRONT 

We suppose that the wave is descending (fig. 10). Denoting the 
slope of the front by 0, j3=— -tan <p. Since we are only considering 
waves whose profiles lie above the free surface, |3 is always positive. 
Denoting, at #=0, the position of the front of the wave by x , and the 
value of the slope by ft>, the equation of the wave profile is thus, 

h=p (x — x), x^x 0> t=0; 



h=0, 


X^Xq, t = 0. 


Equation 112 then becomes 




po(x —x) 2 
°~ 2 ' 


xf£x 0) t=0; 


(7 = 


x^x , t=0. 



(130) 



Introducing the equation for the velocity of propagation of a wave 
volume element, eq 128, we obtain from the above equations 



CO 

If we write 



= ^ 1+ (^ ft "fe) (ab " a!) } *- Xo ' t=0 - 



«b=^- I «oA>-4' (131) 



then the last equation becomes 

co— u =6 (x o — x), x^X ,t=0. (132) 
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Figure 10. — Deformation of the linear front of a descending wave. 

Now let the position of the front of the wave at i=At be x ' (fig. 10). 
The section which at time £=0 had the volume a- in front of it has 
moved from a; to x', and the height has changed from h to h'. Since 
the wave front moves with the velocity co (section III-4), while the 
section at x moves with velocity w, we have 

Xo f —x'=(a) Q —a))At J r(x —x)j x<^x , t = 0. 

Using eq 132, we obtain 

Xo'—x'=(l — d Af)(xQ—x), x<x , t=0. (133) 

But the volume <t in front of x is equal to the volume a t in front of x', 
and thus we have from eq 130 and 133 

<T 1 =^(x '-x'y(l+2d At). 

Comparison with eq 130 shows that during the time A* the linear 
character of the wave front is preserved, but the slope has changed 
from p Q to 0', 



Putting 



we have from eq 134 



j8'=(l+20 o AOA>. 



*'=*+>, 



^=20 O A>, J = 0, 



(134) 



(135) 



which enables us to determine the variation of the slope. Since /3 
is positive, we find that if O >O, the wave is becoming steeper, if 
flo^O, the slope remains constant, if d Q <^0, the wave is flattening out. 
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It is desirable to translate these conditions into relations involving 
P Q . i, and A. Using eq 106, 126, 62, and 48, eq 131 becomes 

fc-f^&fo-D, (136) 

where I—i ( "o-v^— q )* Since g, H , and p are positive, we thus find 

that if &>/, 0o>O; if p =I, 6 =0; and if A><7, <0. Accordingly, 
we rind that, 

0o]>( «•-%/ X/2t— q », the wave is becoming steeper; 
p Q =(-*J\/2i—-]i, the slope remains constant; 

$><(«>-%/ V^i— q V, the wave is flattening out. 

We now derive an expression for /3 as function of time. Since it 
has been demonstrated that the wave front conserves its linear char- 
acter, eq 135 holds for any value of t, and we obtain, using eq 135 
and 136 and discarding the subscript which refers to 2 = 0, 

dp 

dt r 

The solution of this differential equation for the initial condition 
P=p at *=0 is 

(k-ifk-D exploit (138) 

If 1=0, the solution of eq 137 becomes 



l= iylso {0> ~ Ip) ' (137) 



l ~Uk M (139) 



and the slope continually increases. 

We next suppose that the wave is ascending (fig. 11). Denoting 
the slope of the front by p, P=tem <p. We again suppose that the 
wave profile lies above the free surface, which makes /3 positive. 
Using the same method of reasoning as in the above, we obtain cor- 
responding to eq 135 and 136, 



dp 
dt 



=2d Q p { 



0, 



£{-J+M, 



where J 






Translation Waves in Open Channels 
Thus, if ft/lW^ + Qk the wave is becoming steeper; 

^° V 3 A/ 2l^"Q r> ^ ie s ^°P e r ^ m ains constant; 
^(W^'^Ir ^ ie wave * s flattening out. 
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Figure 11. — Deformation of the linear front of an ascending wave. 
The expression for /3 as function of time, becomes 



0— 



A>— (ft— J)exp 



Wfr 7 ' 



(140) 



and if J=0, which for ascending waves can only be the case if the 
slope, i, vanishes, the latter reduces again to 



& 



~2-\H^ 



1 



(141) 



and the slope of the wave is continually increasing. 

If it is desired to express /3 as a function of the distance x ' — x tra- 
versed by the wave, this is obtained by introducing (x '— x )/cc in 
place of t in the various expressions giving fi as a function of t, since 
the head of wave is traveling with the velocity co . 

7. CHANGE OF HEIGHT OF AN ABRUPT WAVE FRONT 

In a wave having an abrupt wave front there is a discontinuity of 
height at the wave front (fig. 12). If for some considerable distance 
behind the front, the height is constant, the change of height with 
time, or preferably with the displacement of the front, can be readily 
determined, using the Boussinesq formula for the velocity of prop- 
agation of a volume element. 
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Fi6 .kk 12,-De/or— 0/ an «*** — *- 

rf that at t-O, in the initial portion of the wave, 

fc=0, x^xo, <=°- 
Since the wave is descending^ 112 becomes 

.•' «f « wave-volume element is thus from 
The velocity of propagation of a wave vo 
eq 128, 106, 62 

From to « •*- « *• ^ the we tro "'' * " (143) 

Wl =coo^H- 4 Wo Ho J 



have 

and from eq 113 






$--»& •<*'-°' 



(145) 



Incomes (I 46 ) 



CM 
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We consider two distinct values of x at 2=0, both behind the wave 
front, denoting them by X\ and x 2) and write hi=h(xi), h 2 =h{x 2 ). 
After an interval of time A/, let the height of the corresponding wave 
elements be W and h 2 . Now 

V— h 1 =-nAt=— t fih Q At, 



and 

dt 



h 2 '—h 2 =-rAt= —fihoAt. 



Since hi=h 2 =h Q , we see that hi'=h 2 '. This shows that the con- 
stancy of height of the wave front is preserved. Replacing h Q by h in 
eq 143, we see that the velocity of the wave front a>i is 

Again replacing h by h, eq 146 now becomes 

j t =—fih, x<x + I o)!^. (148) 






Integrating this equation, and using the initial condition h=ho at 
2=0, we obtain h as a function of t 



h=h Q e ht , x<^x -\- o)idt. 



(149) 



If it is desired to determine h as a function of the distance Xq' — x q 
traveled by the wave front, we write from eq 147 and 148 

®L =m x -l^M£l (150) 

dx w L 4 w HqJ 

since dx 'ldt=oo 1 ', and 3h^gH /4:co H is required to be small with 
respect to unity. The solution of this differential equation for the 
initial condition h=h j at x Q / —x is 

3 -y/gHo h 

4 co i? 

8. SOLUTION OF THE WAVE EQUATION WHEN SLOPE AND FRICTION 
ARE NOT NEGLIGIBLE 

The significance of introducing w, the celerity of a wave-volume 
element, lies in the fact that the integration of the wave equation, eq 
111 becomes equivalent to solving 

g+!<W-ft (113) 
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where co is given by eq 116, or in this particular case by eq 128. 
Although the general solution is complicated, if it be assumed that 
b 2 h/dx 2 is negligible, compared to 6A/i? 2 ^2, the problem is consider- 
ably simplified. This assumption requires that the curvature of the 
wave profile be small. We thus obtain from eq 129 and eq 62: 

uh=hU + (co - U*)(l +^ )^=F rfu (152) 

and eq 113 becomes 

dh 



s+" 



-[*7 + (coo- ^o)(l +^ )]§= ~¥i, (153) 

since do-/dx= T h. The general solution of this differential equation is 

x=(u + («o- U )^ ^l-iy+sihe'V, (154) 

where <£ is an arbitrary function. This integral is due to Boussinesq 
[3, p. 457, eq 392] and Favre speaks of it as one of the Boussinesq 
integrals [6, p. 106]. 

In determining the arbitrary function <f>, two cases are of particular 
interest; (1) the channel is infinite in both directions, and (2) the 
channel has an entrance and is infinite in only one direction. In the 
first case, if for t=0, h is a known function of x, — o><^x<^cq } then $ 
can be obtained directly. Let the function be h=^f(x) f t=0. Put- 
ting £=0 in eq 154, we have x— $(A), #=0, hence in this case <i> is 
simply the inverse of ^. In the second case, $ is determined from 
the boundary condition &=#(*) at the channel entrance x=Q, but in 
this case the evaluation is indirect and tedious. At #=0, eq 154 
gives 

*(a«v) = -["o+ (coo- ^mr^r^Y (155 > 

In general, the form of the equation does not permit solving explicitly 
for either t or A, hence it is necessary to evaluate $ numerically. By 
considering values of t in the range — °° to + °o , and the corresponding 
values of h at the entrance, using h=i[/(t), the magnitude of the 
independent variable can be computed. Corresponding values of $ 
are obtained from the right-hand member of eq 155 by inserting the 
same values of h and t used to evaluate he tf i. 

9. INTUMESCENCES PRODUCED BY OPERATING LOCKS IN A CANAL 

CARRYING WATER 

We now consider a practical example of case (2) above which can be 
treated by a simpler method. In many cases the canals leading to 
and from hydroelectric plants are also used for navigation, and a 
lock is provided to bridge the difference in head utilized by the plant. 
Whenever this lock is emptied, a positive wave is produced in the 
downstream canal, and, whenever it is filled, a negative wave is pro- 
duced in the upstream canal. The former, of course, travels down- 
stream, whereas the latter travels upstream. 
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Because of such design considerations as the possibility of over- 
topping the embankments, it is desirable to predict by computation 
the magnitude and deformation of these waves. Ordinarily the 
discharge from the lock is a continuous function of time, the exact 
relation being obtained from computation or model experiments. The 
waves produced by operating the lock are small in height and long 
enough so that the curvature is negligible, and thus eq 154 is applicable. 

As far as this analysis is concerned, the lock constitutes a termina- 
tion of the canal, and the boundary conditions of case (2) apply. 
The complications of this case are eliminated by a method due to 
Favre [6, p. 112], which is based on the fact that the distance traveled 
by the wave is very large compared to its length. The initial shape of 
the wave is computed for the initial time #=0, in the form h=^(x) } 
assuming that the effect of friction is negligible during the short time 
interval required for the formation of the wave. This function is then 
used to determine <£. 

Let the discharge at the lock per unit width of channel be represented 
by the function 

2-/(0, x=0. (156) 

As a result of the above, we may also write 

f**#(2), x=0. (157) 

It follows from eq 154 that the form of the wave at the moment 
it is produced and shortly thereafter, neglecting the effect of friction 
during the genesis of the wave, is represented by 

X=[co +(co -?7o)^J + ^W, (158) 

since Lim — *^z? -=1. From eq 52 and 117 the mean velocity 

tfi->0 VI 

in a cross section beneath the wave is 

u=Uo+ wn h - (159) 

Neglecting the effects of friction and curvature in eq 129, this becomes, 
using eq 62 

u = u ^^- u K 1+ m)H^h- (160) 

The discharge per unit width at x=0 is thus 

q=HTJ= U (H +h) + (o>- U )(l+^h=f(t), x=0. 
Comparing this equation with eq 157 we have 

t=F{Uo(H +h) + (cjo-U )(l+^yi) f 3=0, (161) 

524586—43—6 
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which expresses h as a function of t for £=0, and thus serves to deter- 
mine the arbitrary function $ in eq 158. Placing x=0 in that equa- 
tion, we have 

*(A) = -[**+ (co - Uo)^, < 162 ) 

and from eq 161 

$(/>) = -[«*+ (co»- U^JfiUoiHo+h) + ( w „- U )(l +tf£)V- (163) 

Since the function $ in eq 158 is the same as the function $ in eq 154, 
we write for the second term in the right-hand member of this equation, 



$Qie 



ht ) = -[uo+(o> -Uo) 



3hef> 



hi-] 

+ (o> -U )(l 






(164) 



The solution of a typical problem will now be described in detail. 
It is desired to find the wave profile in the canal at any instant t when 
the discharge is given as a function of time, Q/b=q=f(t) J b being the 
width of the free surface of the canal. We shall assume that the 
discharge first increases and then decreases, as shown in figure 13. 
This will create a positive surge in the discharge canal and a negative 
surge in the supply canal. The origin in figure 13 is conveniently 
chosen to correspond to the instant of maximum discharge. The 
abscissa t is, of course, expressed in seconds, and the ordinate, q 7 in 
any convenient unit of discharge per unit width, such as square feet per 
second. In eq 154 and 164, upon which the solution depends, U and 
co will have to be expressed in corresponding units of velocity, such 
as feet per second. 




Fig. 13. — Discharge from supply canal as a function of time: Q/b=q—f(t) y eq 156 ', 

or t=F{q) y eq 157. 

The positive surge in the discharge canal will now be considered. 
The quantities U , H , R and i are determined from measurements of 
the steady flow. The coefficient of friction, Xp, is then determined from 
eq 48. AH the quantities needed to determine co from eq 62, and j x 
from eq 126 are thus available. Note that the positive sign in eq 62 
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is selected, since the wave will be moving downstream. A small posi- 
tive value of h is now selected, which we shall denote as hi. The argu- 
ment of the function F in eq 164 is computed, using the value of t 
corresponding to the instant for which the wave profile is desired. As 
is seen from eq 157, this gives a corresponding value of g, which we shall 
call ft. If g x is greater than the maximum value of q which is actually 
obtained, this signifies that hi was higher than any point of the intu- 
mescence. It will be advisable to select hi very small at first and then 
increase it until the top of the curve in figure 13 is attained. The values 
of t corresponding to hi will be denoted by t n and t n since, there will 
generally be two of these, the first positive and the second negative. 
This is clearly indicated in figure 13. These two values, tn and t i2 , are 
thus the desired values of the function F occurring in eq 164, and the 
two values of the double-valued function $ can now be obtained by 
replacing h in the coefficient of Fbj the quantity hi. Turning now to 
eq 154, replacing h by hi and t by the value corresponding to the instant 
for which the wave profile is desired, we obtain two values for x } cor- 
responding to the original value h u which we shall denote by x n and 
x 12 . These are the abscissas of the points on the profile having the 
height hi (see fig. 14). The whole process is repeated for greater 
heights, h 2 , h s , . . . , until sufficient points are determined to estab- 
lish the contour of the wave. 

The process for determining the negative surge in the supply chan- 
nel is exactly the same, except that in eq 62 the minus sign is chosen, 
since the wave is traveling upstream. The quantities h u .... and 
#n, #12, . . . , will be negative, as shown in figure 15. 




Figure 14:.— Profile of descending wave in discharge canal. 



10. EFFECT OF FRICTION ON WAVES WHEN THE SLOPE OF THE 

BED VANISHES 

We now consider ease (II) of section IV-5, i=0 and U =0. The 
determination of the positive sense of the z-axis on the basis of the con- 
ventions in section III— 1 fails, since there is no flow. Furthermore, 
the terms "ascending" and "descending" waves explained in section 
III-3 fail to have any significance, since we have no method of dis- 
tinguishing the upstream from the downstream direction. We thus 
are forced to make an arbitrary choice of the positive z-direction, and 
consider downstream to be in this direction, A descending wave thus 
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Figure 15. — Profile of ascending wave in supply canal. 

becomes one for which co >0, and an ascending wave one for which 
co <0, that is to say, a descending wave is propagated in the positive 
^-direction, and an ascending wave is propagated in the opposite 
direction, since our hypothetical observer is stationary and the down- 
stream direction coincides with the arbitrarily chosen positive sense 
of the #-axis. 

Since i=0 and U =0, it follows immediately from eq 52 that 



gHi—^ZP^—^uK 



(165) 



It has previously been tacitly assumed that U was always positive, 
that is, u was smaller than U , so that U and U had the same sign. 
Since in this case U =0, this assumption is no longer valid. The 
sense of the friction is always opposite to that of the velocity, and thus 
if u is negative, the sign of the right-hand member of the above 
equation must be changed, since we wish to keep X positive, regardless 
of the sign of t . From the first approximation, eq 64, u=co h/H Q , and 
thus the velocity must have the same sign as oo h. We introduce the 
new parameter 



J2— 4> 



if oj /i>0; 



/ 2 =-^ ,if^<0; 



and obtain from eq 165 



gHi-^U 2 -- 



:-2co o2 (|)/, 



(166) 



(167) 



The difference between \ u 2 and \u 2 is a quantity of the order being 
neglected in this analysis. If we introduce the quantity 



G '=± \h 2 dx, 



(168) 
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then since da'/dx= Th 2 } we can write eq 167 in the form 

0Hi-lu^±2^(-Q (169) 

Substituting this in eq 116, we have for the velocity of propagation 
of a volume element 



oo = ai ( 14 







2+ki h 



H a 



k 2 H 2 d 2 h. 
Qh dx 2 ' 



w — U hH 2 



«0 



> 



(170) 



Since in this case Z7 — 0, we see from eq 106 that ki=k 2 =l, and thus 






6\ 

oV 



(171) 



Consider the sign of o) a[f 2 /hH 2 y which equals (a^&'fj/onJiHo 2 . Since 
from eq 168, <r' is always positive, it follows that the sign of the above 
is the same as that of / 2 /co A. But from eq 166 f 2 has the same sign as 
wji, hence the original term is always positive. In accordance with 
the convention of choice of signs for descending and ascending waves, 
the choice in eq 171 is always opposite to that of the sign of a> . There- 
fore, the effect of friction on the celerity co when i=0, U Q = 0, is to 
diminish its absolute value by an amount which increases with the 
distance from the head of the wave. 



11. SOLUTION OF THE WAVE EQUATION WHEN THE SLOPE 
VANISHES BUT FRICTION IS CONSIDERED 



As in section IV-3, the problem is to solve 



(113) 



where u is given by eq 116 or, in this particular case, by eq 171. For 
the reasons mentioned in the same section, we assume that b 2 h/dx^ is 
negligible. Equations 171 and 113 thus become 

(172) 



and 



/.. 3h _ <r'/ 2 \ 

bh . /,-, U\dh . /h V 

Tt +Uo { 1+ 2H hx = -J 2W \H ) ■ 



(173) 



respectively, since da'/dx=Th 2 . Assuming that the wave is every- 
where positive or everywhere negative, then / 2 is constant, and the 
general integral of the above differential equation is 



£ = COo 



1+ 



^iog(i-^y 



2H 



'Ho 2 



TT 2"^ 



r h i 



(174) 



where $ is an arbitrary function. This is Boussinesq's second integral 
[3, p. 454 eq 393]. 

Regarding the determination of this arbitrary function, we consider 
the same two cases as in section IV-8: (1) the channel is infinite in 
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both directions, and (2) the channel has an entrance which is infinite 
in only one direction. In the first case, if for t=0, h is a known func- 
tion of x, — oo <z< oo 9 then <£ can be obtained directly. Let the func- 
tion be h=y(x), t=0. Putting t=0 in eq 174 we have x=$(h), t=0, 
hence in this case $ is simply the inverse of ^. In the second case, <f> 
is determined from the boundary condition h=\l/(t) at the channel 
entrance x=0; but in this case the evaluation is indirect and tedious. 
At x=0, eq 174 gives 



<£/ -? 1 = — w 



V3w 



3h 



i<, g (i-^y 



2H — /2W0 1 ■ . 
H 2 m 



(175) 



In general, the form of the equation does not permit solving explicitly 
for either t or h, hence it is necessary to evaluate <2> numerically. By 
considering values of t in the range — <» to +» , and the corresponding 
value of h at the entrance, using h=\l/(t), the magnitude of the inde- 
pendent variable can be computed. Corresponding values of $ are 
obtained from the right-hand member of eq 175 by inserting the same 
values of h and t used to evaluate 

h 



1 H jit 



These methods are precisely analogous to those discussed in connection 
with eq 155. 

12. INTUMESCENCE PRODUCED BY OPERATING LOCKS IN AN 

ORDINARY CANAL 

In contrast to canalized rivers and hydroelectric canals, such as 
discussed in section IV-9, ordinary canals contain still water except 
when disturbances are set up by the operation of the locks. The 
general considerations at the beginning of section IV-9 (through eq 
157) also apply to this case, except that i=0 and t7o=0, and hence 
we must consider case (2) of eq 174. 

The discharge per unit width at the lock is again represented by 

2=/W, x=0, m (156) 

or by the inverse function 

t=F(q), s=0. (157) 

Neglecting the effect of friction during the formation of the wave 
and shortly thereafter, we obtain from eq 174 

a=o*(l+^ o )*+*(A), (176) 



since Lim -? — = 1 . 



los(lJp*<) 
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For the mean velocity in a cross section beneath the wave, eq 52, 
117, 171 give 

o^+Sdnhv- (177) 

since U =0, and friction and curvature are neglected. The discharge 
per unit width at x=0 is, therefore, 

(l=HU=uJi(l+~^\ 3=0. (178) 

Substituting this value of q in eq 157, we obtain 

t=F(a h(l +||;))> x-0. (179) 

But from eq 176 

*(A) = -«o(l+^)*, *=0, (180) 

or from eq 179 

•w— ^i+^X**( 1+ S)} (181) 

Since this function <£ in eq 176 is the same as the function $ in eq 
174, we have for the right-hand member of this equation, 

/ * \ r 1+ 5* 

i^-sA H iJ -sm <182) 

Equations 174 and 182 are used to determine the wave profile at 
any instant t, as is described in section IV-9. 

13. LIMITATIONS OF BOUSSINESQ SOLUTIONS 

It is necessary to emphasize the fact that the second-order theory, 
which is based on the concept of the velocity of propagation, co, of a 
volume element, is subject to certain important restrictions. The 
most general expression for o> is given in eq 116. In deriving this 
equation, the important assumption was made that oo/o>o is nearly 
equal to unity. 

We shall now turn our attention to the specialized equations derived 
from eq 116. These are eq 118, 128, 171. Examination shows that 
in each case, w/co =l + (correction terms). The restriction on the 
theory thus requires that these terms shall be small compared to 
unity. The first term, which is common to all equations, is (2+ki)h 
/4H . This merely requires that h^4H /(2+ki). The second term, 
which expresses the effect of curvature, requires that (b 2 h/c>x 2 )/h 
<^.6/H 2 k 2 . We see from this, that (b 2 h/(>x 2 ) must vanish to a higher 
order than h. 

The third term, which appears only in eq 128 and 171, expresses the 
effect of friction, These two equations are applicable only when 



$ 
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aji/o) h or cr'f 2 lhH 2 , respectively, is small. But a and a' can be made 
as large as we please merely by considering longer and longer intumes- 
cences, and if the intumescence has a finite volume, h approaches zero 
at the tail whereas a and cr / remain finite. This shows that the term 
under consideration may be large, or even infinite. In view of this 
fact, it may be expected that these solutions will give correct values 
only in the frontal portion of the wave where <s\h and cr'/h are small. 
It should be noted that this remark applies with equal force to the 
Boussinesq integrals, eq 154 and 174. 
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